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1. Introduction 



In Quantum Chromodynamics (QCD), the ground state must be the vacuum in order to 
ensure the cluster decomposition property of physical observables. However, in unquenched lattice 
QCD simulations by Hybrid Monte Carlo algorithm [|lj], the correct sampling of topological charge 
will become increasingly more difficult In view of this situation, one promising approach is 

to fix the topology during the Hybrid Monte Carlo simulation and try to extract physics from the 
simulations at a fixed topological charge Q. It is known that the fixed Q effect is a finite size effect 
which vanishes in the infinite volume limit. In this report, we give a theoretical basis for estimating 
the finite size effect in order to extract physics in the vacuum from QCD at fixed topology [Q] , 
by extending the work by Brower et al. We derive a general formula to estimate this relation. 
Furthermore, using our formula, we propose a method to measure the topological susceptibility at 
fixed topology. The numerical study is given in a separate reports [§, 0]. 



2. General formula 



Consider the partition function in the vacuum defined by 

z(e) = (e\e) = exp[-vE(e)], (2.1) 

where E(6) is the vacuum energy density. The topological susceptibility % t at = is defined by 

(0|<2 2 |0) d 2 E(d 



Xt 



V dd 2 



(2.2) 

=o 



Since % t > by definition, = is a local minimum of E(6). Moreover, Vafa and Witten proved 
that Z(0) > Z(9) [|8[], so that = is the global minimum of the function E(6). Assuming 
analyticity of E(6) near = 0, we can expand E(6) as 

= L 7^y d2n = f 2 + O(d 4 ). (2.3) 
The partition function at a fixed topological charge Q is a Fourier transformation of Z(6) 

Z 2 = T- [*d0Z(0)cxp(i0Q) = — [* ddexp(-VF(d)), (2.4) 

2K J -71 2K J -71 



where F(6) = E(6) — idQ/V . For a large enough volume, we can evaluate the integral in (2.4) 
by the saddle point expansion. The saddle point 6 C is given by 

e c = i^(l + 0(8 2 )), (2.5) 

where 8 = Q/ (XtV)- We then expand F(6) as 

E^i „ ^ EW(9 ) 

F(d) = F(d c ) + — (0 C )(0 - C ) 2 + £ ^(0 - C )», (2.6) 

2 , n\ 
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where Z?W is the n-th derivative of E(6) with respect to 6 at 6 = 6 C , and is given by 



VF(d c ) = + 0(8 2 )), E( 2 \e c ) = X r(l + 0(5 2 )), 

E {2n \d c ) = c 2n (l + 0(S 2 )), E^ 1 )^) = 6 c c 2n (l + 0(8 2 )). 
By a change of variable s = VE^V (6 — C ) we can rewrite the integral as 

,2 E {n)y 



VF{B C ) r VWfv{iz-e c ) 



ds exp 



2 „= 3 " ! VV^V 



(2.7) 
(2.8) 

(2.9) 



Neglecting exponentially suppressed terms and expanding in powers of 1/V, we obtain 



1 



exp 



Q 2 



2XtV 



(2.10) 



This shows that, as long as 8 <ti 1 (equivalently, Q ^ &V), the distribution of 2 becomes the 
Gaussian distribution. Similarly, consider an arbitrary correlation function in 6 vacuum and at 
fixed topological charge Q are defined as 



G{d) = (d\Oi0 2 ---o n \o), g q 



1 1 

Zq 2k 



ddZ(d)G(6)exp(idQ). 



(2.11) 



Using the saddle point expansion as before, if G is CP-even, we can show 

Q 2 



G ^ = G (Q )+G ^ {0) _L_ 



1 



X,V 2 X ?V 



+ c (4) (0)^ + o(v- 3 ) 



(2.12) 



where G" (0) stands for the n-th derivative of G with respect to 6 . If G is CP-odd, we have 

Gt = G(1)(0) J_ (l-^) +G (3) (0 )-^ + O( y-3). (2.13) 



The formula ( [2.1 2[ ) provides an estimate of the finite size effect due to the fixed topological 
charge. The leading correction is of order 0(1 /V). It should be also noted that the other formula 
( 2.13 ) suggests that it is possible to extract the 6 dependence of CP-odd observables, such as the 
neutron electric dipole moment by measuring the observable at a fixed non-zero topological charge, 
once the topological susceptibility % t is obtained. 



3. Topological susceptibility 

There are two reasons for measuring topological susceptibility. The primary reason is to study 
whether the local topological fluctuation is sufficiently created in unquenched QCD simulations. 
In conventional quenched QCD simulations, the topology change during Monte Carlo updates 
is triggered by the formation of dislocations which grows into local topological fluctuation with 
positive or negative topological charge. This is a topology non-conserving processes through lattice 
artifact. On the other hand, in unquenched simulations such process is highly suppressed towards 
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continuum and chiral limit. In particular, unquenched simulations by the JLQCD collaboration 
explicitly prohibits such processes by introducing extra Wilson fermions Then, the main source 
of the local topological fluctuation is achieved through the pair creation of local fluctuation of 
positive and negative topological charge densities, just as in the continuum theory. Therefore, 
the measurement of the topological susceptibility is a crucial test of the thermal equilibrium in 
local topological fluctuation. The secondary reason is that the topological susceptibility is the key 
quantity to estimate the finite size effects at fixed topology as was shown in the previous section. 

Recently, there are several proposals for the field theoretical definitions of the topological 
susceptibility, which are free from ambiguities. The first one is given by Giusti et al. [ |10| ] in 
which they define the topological susceptibility by the integration of disconnected flavor-singlet 



pseudoscalar correlator. The second one is a UV divergence free definition by Luscher [11] in 



terms of n-point function of flavor non-singlet scalar and pseudoscalar correlator. The third one is 



proposed in the study of Schwinger model [12]. They extracted the topological susceptibility from 



the asymptotical value of the singlet pseudoscalar correlator up to 1 /V correction as 

lim (mP(x)mP(0)} Qy = ±-(Q-- Xt ) + ^(y- 2 ). (3.1) 
x->°° V \ V ) 



Although the intuitive picture for this relation was given in Ref. [|12[], the field theoretical proof 



based was given only recently, which will be explained in the next subsection. 

3.1 Field theoretical proof of the formula for the topological susceptibility 

Suppose that there is a well-defined local operator co(x) that measures the local topological 
charge. The global topological charge Q is then obtained as Q = / d 4 xco(x), and the topological 
susceptibility is % t = / d 4 x((o(x)co(0)), where the expectation value is taken for the 6 = vacuum. 
Since (o(x)a>(0) is CP-even, Eq. ( |2.12[ ) gives 

1 A Q 2 c 4 



<fl)(jc)fl)(0)) e = <o(jc)fl)(0)> + (co(x)co(0))W ^— i 



2V Xt V VXt 2 Z ?V 



1 



+( (O (x)(o(0))^— — + 0(V- i ), (3.2) 

where (O)W is the n-th derivative of (O) with respect to 6. In the large separation limit \x\ — > oo, 
the CP invariance at 6 = and the clustering property at a fixed 6 gives 

hm (o)(x)o)(0)) G = i f f - Xt -^-) +0 (y- 3 ) + 0(e-^), (3.3) 

where the flavor singlet pseudo-scalar meson mass, m n >, is the lightest mass of possible intermedi- 
ate states. 



Physical quantities such as the topological susceptibility % t can be obtained through (pT3|). In 
practice, this formula will be used for a finite separation x instead of \x\ — > °°. The clustering 
property in the 6 vacuum receives a correction of order of g _m i'W, which vanishes quickly because 
the flavor singlet meson r\' acquires a large mass due to the axial anomaly of QCD. 

We now express the bosonic correlation function (co(x)co(0)) in terms of a fermionic one using 
the anomalous axial U(l) Ward-Takahashi (WT) identities for an arbitrary operator O: 

{d ll A tl (x)0-2mP(x)0 + 2co(x)0 + 8 x O) = 0, (3.4) 
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Figure 1: The finite size effects 1 + T m (left panel) and 1+77 (right panel) from fixed topology at Q = with 
L = 2 fm. The pion masses correspond to those in the dynamical simulation by the JLQCD collaboration. 



where (x) 



W~ £/ ( x )7^Y5¥ f (x) and P(x) = ^'£fW f (x)YsY f M are the flavor singlet axial- 
vector current and pseudo-scalar density, respectively, and 8 x O denotes a axial rotation of the oper- 
ator O at x. Combining WT identities for O = 2mP(Q) and O = 2(0 (x) Combining these we finally 
obtain 

l ( Q 2 



Urn (mP(x)mP(p)) Q = ~(2--Xt 

l^large V \ V 



C4 

2X,V 



+ 0(e 



(3.5) 



In fact, using this formula one can determine the topological susceptibility from the nonzero asymp- 
totic value of the singlet pseudoscalar correlator [§, ^]. 

4. Application to physics 

One can estimate the finite size corrections for the pionic quantities with the help of Chiral 
Perturbation Theory (ChPT). Using the next-to-leading order ChPT formula, 6 dependence of the 
the pion mass and the decay constant for two-flavor QCD is 



ml(d) 



Me) 



\rmNLO 



\rmNLO 



ml(9) 



f 




(4.1) 
(4.2) 



where m n (Q) 



2Bom q cos an d ^3, U are the low energy constants which can be estimated as 

ff ys = 2.9±2.4, and/J hys =4.4±0.2 [|l7|]. Fig. [T] shows the finite size effects 1 + T m = mf°/% 
and 1 + T f = f^°/fn at L = 2 fm using the NLO ChPT. The pion masses correspond to those in 



the dynamical simulation by the JLQCD collaboration [18]. The topological susceptibility % t is 
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extracted from the singlet pseudoscalar correlator as explained in the previous section. It is found 
that finite size effect for the pion mass ranges from 0.5 to 2.5 %, whereas that for the pion decay 
constant is well below 0.5%. These finite size correction from fixing the topology is correctly 
taken into account in the spectrum study [18]. In general, for quantities which has a non- vanishing 
chiral limit, 6 or Q dependence correction only comes through m n (0) as sub-leading corrections. 
Therefore, pion receives the largest finite size correction. This means that if the finite size effect of 
the pion mass in under control, other hadronic quantities are safe. 



4.1 Nongaussianity of the topological charge distribution 



Recently deviation of the topological charge distribution from Gaussian is observed for pure 
Yang -Mills gauge theory [|l3|, [l4|, 15j 16]. In principle, we can also measure the deviation from 
the Gaussian (C4 coefficient) by combining the 2-, 3-, 4-point functions of the topological charge 
density given as follows 



lim (fl>(x)ffl(0)) e 

br — »large 



lim {(o(x\)(o(x 2 )(o{x'i))Q 
\xj— Xj\— >large 

lim {(o(x\)(o(x2)(o{x^)w(xa))q 

|*r-jt/|-vlaige 



Xt_ 

V 2 



1 



2X?V 



(c 4 -2XtQ 2 ) 



+ o(v- 



(4.3) 



-3*, 

T Xt 

V 



Q 

*V2 



1 + 
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(c 4 --XtQ 2 



+ 0(V- 4 ) (4.4) 



+ 0(V 



-4\ 



(4.5) 



5. Summary 

We have derived general formulas which express arbitrary correlation functions at a fixed 
topological charge Q in terms of the same correlation function (and its derivatives) in the 6 vacuum. 
The difference between the fixed Q vacuum and the fixed 6 vacuum can be shown to disappear in 
the large volume limit as 1 /V only using fundamental properties of the quantum field theory. 

These formulas open a new possibility to calculate physical quantities in the lattice QCD 
simulations at a fixed topological charge. This will become unavoidable as the continuum limit is 
approached, irrespective of the lattice fermion formulation one employs, as far as the algorithm is 
based on the continuous evolution of the gauge field. 

Applying the formula for rc-point functions of the topological charge density ft), we have shown 
that the topological susceptibility %, and C4 appear at the first and second order corrections, respec- 
tively. In principle, these parameters can be determined by the lattice data. Our method is free from 
short-distance singularities, since the local topological charge operators are put apart from others 
and no contact term appears. Numerical calculation is in progress by the JLQCD collaboration 



on the gauge configurations generated with dynamical overlap fermion [fljj, 20, 21, E2p. Once 



these parameters are numerically obtained, they can be used as input parameters for other physical 



observables, such as weak matrix elements such as Bk [23], pion form factor [[24]], the neutron 
electric dipole moment, and so on. 
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